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ABSTRACT 



A storage medium encoded with machine-readable com- 
puter program is used for finding the shortest path between 
two points on a polygonal surface. The storage medium 
includes instructions for causing a computer to implement a 
method for finding the shortest path. Given a first point and 
a second point on the polygonal surface, a polyline lying on 
the surface and passing through the two points is defined. 
The polyline on a polygonal mesh is analyzed to determine 
points lying on the polyline and on edges of the mesh. The 
polygonal faces of the mesh are assumed, without any loss 
of generality, to be triangles. If the start and end points of the 
polyline are not on the edges of the mesh, the faces of the 
polygonal surface on which the start and end points of the 
polyline lie are triangulated so that the start and end points 
become vertices of the polygonal mesh. The polyline is then 
modified such that it will pass through the first and second 
points of the polyline, creating a new polyline of a shorter 
length. The analysis, triangulation and modification process 
are repeated iteratively until a shortest possible polyline is 
found between the first and second points. 

20 Claims, 3 Drawing Sheets 
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ITERATIVE DETERMINATION OF THE 
SHORTEST PATH BETWEEN TWO POINTS 
ON A POLYGONAL SURFACE 

BACKGROUND OF THE INVENTION 5 

This invention relates to a solution for the local minimum 
length path drawn on the surface of a polyhedron opened or 
closed surface. 

Three-dimensional surface renderings have commonly 
been made using a constant triangular mesh over the surface 
of interest. That is, a surface of interest of a computer- 
generated three-dimensional image representation is 
selected. For example, one might use a nuclear, CT 
(computerized tomography) scanner, or other diagnostic 
imaging device to generate a three-dimensional image rep- 
resentation of a portion of a patient's torso surrounding the 
liver. Using conventional techniques, such as comparing the 
gray scale or CT number of each pixel in the image with a 
characteristic gray scale or CT number for the liver, one can 2Q 
identify the liver and the surface of the liver. Because all of 
the liver surface tissue has substantially the same gray scale 
or CT number, if a perspective view of the liver were 
re-created and displayed on a video monitor, the resultant 
display would have the appearance of a silhouette. There 25 
would be substantially no surface texture information con- 
veyed. 

In order to convey contour information, the surface of the 
liver or other structure which is visible from a viewing 
direction is overlaid with a constant triangular mesh. That is, 30 
the surface of the liver or other object of interest is approxi- 
mated by a mesh of triangles of as uniform a size as 
permitted by the surface texture. In a typical 512 times 512 
frame for a medical diagnostic image, one might expect to 
find about 20,000-50,000 triangles on the surface of interest, 35 
The vector surface normals are each compared with an 
illumination direction vector and a viewing direction vector. 
Based on the vector comparison, a gray scale (or color or 
hue) for each triangle is selected. Typically, the gray scales 
are selected to approximate the lighting and shading condi- 40 
tions that an object of the selected configuration would have 
at the corresponding point on the surface when viewed from 
the viewing direction and illuminated from the illumination 
direction. 

The display may represent more than just the surface of 45 
the object. In particular, one can commonly define a cut 
plane on a medical diagnostic image and remove the portion 
of the object to one side of the cut plane. The portions of the 
image corresponding to the interior of the liver or other 
object of interest are displayed with appropriate gray scales 50 
or colors which correspond to the CT numbers or pixel 
values of each internal pixel of the liver or other structure 
that is intercepted by the slice. The remaining surface 
portion continues to be displayed with three-dimensional 
surface rendering. 55 

Typically, the operator has the ability to rotate the dis- 
played object. Each time the object rotates an incremental 
distance, the surface normals are again computed and com- 
pared with the illumination and viewing direction vectors 
and reassigned gray scale values to each triangle in accor- 60 
dance with the comparison. The large number of calcula- 
tions necessary to reassign the gray scales are commonly 
performed slowly relative to the rotation rate of the object. 
That is, the operator rotates the object faster than the new 
gray scales can be computed, causing the object to appear to 65 
rotate in jumps, rather than smoothly. For example, if an 
object rotates about 5 degrees in the time necessary to 
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recalculate the gray scale values, then the object appears to 
rotate in 5 -degree steps. 

"Decimation of Triangle Meshes" by Schroeder, et al., 
Computer Graphics, Vol. 26, No. 22, pgs. 65-70, July 1992, 
describes an algorithm which reduces the number of tri- 
angles required to model a physical or abstract object. In the 
Schroeder technique, multiple passes are made over an 
existing triangular mesh, using local geometry and topology 
to remove vertices that pass a distance or angle criterion. 
Holes left by the vertex removal are patched using a local 
triangulation process. One difficulty with the Schroeder 
technique is that it is relatively slow. Running time for the 
technique is on the order of minutes. Commonly, radiolo- 
gists expect to review CT, nuclear camera and magnetic 
resonance images substantially in real time and do not find 
waiting on the order of minutes acceptable. Arata U.S. Pat. 
No. 5,689,577, issued Nov. 18, 1997, teaches an improved 
three-dimensional triangle mesh simplification technique 
which overcomes the above referenced problems. 

A common problem in computational science is the need 
to represent the geometry of an object or the structure of 
data. One method for representing such objects is to use a 
polygonal mesh. A system and method that generates a large 
polygonal mesh is described in Cline et al. U.S. Pat. No. 
4,710,876, issued Dec. 1, 1987, and assigned to the instant 
assignee. Other techniques, such as automatic mesh genera- 
tion or terrain mapping from satellite data, are also capable 
of generating large polygonal meshes. For example, the 
complex, curved surface of a human tooth can be approxi- 
mated by using many thousands of triangles (or other 
polygon types) joined along their common edges. The ability 
to represent geometry is important for many reasons. In 
computer graphics, polygonal meshes are used in the light- 
ing and shading operations to generate images. Polygonal 
meshes are used in numerical analyses to represent the 
boundary of solid objects. From these representations, equa- 
tions can be developed to solve such complex problems as 
heat flow, structural response, or electromagnetic propaga- 
tion. Another application is in geometric modeling, where 
polygonal meshes are often used to determine object mass, 
volume, surface area, center of gravity, and moments of 
inertia. 

In the past, polygonal meshes were typically comprised of 
hundreds to thousands of polygons, and computer hardware 
and software has been designed to process the volumes of 
information obtainable therefrom. However, recent 
advances in computational science have resulted in tech- 
niques that generate hundreds of thousands or even millions 
of polygons. Such large numbers, while capturing the geom- 
etry of the object very precisely, often overwhelm computer 
systems. For example, most graphics systems presently are 
incapable of rendering a million polygons at a speed that is 
not detrimental to interactive computation. 

The basic problem is that techniques that generate large 
polygonal meshes are extremely valuable and cannot be 
easily modified to produce fewer polygons. Hence a general 
technique for reducing, or decimating, a mesh composed of 
a large number of polygons to one containing fewer poly- 
gons is necessary. Furthermore, the decimation process, to 
be effective, must preserve the topological and shape prop- 
erties of the original polygonal mesh. 

Zarge et al. U.S. Pat. No. 5,590,248 (the '248 patent), 
issued Dec. 31, 1996 and assigned to the instant assignee, 
teaches a method for reducing the complexity of a polygonal 
mesh while preserving the topological properties and shape 
of the original mesh. By this method, each vertex in the 
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mesh is analyzed to determine if it is superfluous and is erate medical information and display images at the option 

removed along with all of the triangles connected to it if it of the viewer. The viewer is provided with the option of 

is superfluous. The polygon created by removing a vertex is selecting one or more threshold values which determine, for 

filled (retriangulated) with new triangles according to an example, whether or not bone surfaces as opposed to brain 

algorithm disclosed herein. The method further includes the s surface tissue is to be displayed. The viewer or operator of 

step of removing triangles (consolidation) from the mesh if the P reserjt system can also select the appropriate viewing 

the triangles are relatively small compared to neighboring a , n S le and > can at Wll j» selectively ignore segments of the 

triangles. The method also includes removal of edges from d * ta generated m order to provide cross sectional views 

the mesh if those edges are relatively short as compared with trough any desired plane. 

neighboring edges. 10 Because the viewing angle is selectable, it is possible to> 

„_ 4 . . . , , 4 generate a sequence of images and display them sequentially V. 

2f P at f enl V^ Ta , V , C °7 U -£ r Pr °" » P^vide the medical practitioner with interior views of 7 

grammed to perform the method steps disclosed. Tte pro- fa surfaces £ dimensional manner from any-' 

gram allows a user to specify various thresholds therein , . , . . , ... - ..... ru . 

°. . , i i r *u desired viewing angle, with the further capability of being 

disclosed to control the representational accuracy of the . mM# * „, »k™„u „i™ ^ cl ,. M u/l-iT 

. ™ ^ . , • k able to-construct a view through any plane or slice. While 

resulting mesh. The program may also be run in an iterative 15 - . * • a i 

r , . 5 i • , . . l i- iL for many purposes, an almost infinite variety of meaningful 

fashion until a desired reduction in mesh complexity has . ' F A frnm . o ' nf XjtD _ y-r 

. . . . images may be created from only a single set oi MR or CT1 
neen achieved 

scan slice data array, if the objective of the medical inves- 

The aforementioned U.S. Pat. No. 4,710,876 (the '876 tigatioQ is tne study of internal anatomic variations as a 

patent) is generally directed to a system and method for function of time, then it is meaningful to produce a sequence 

displaying surface information. The images of the surfaces of three dimensional data arrays indexed by time. The 

displayed are typically contained within the interior regions system and me thod of the present invention provide the 

of solid bodies which are examined by computed tomogra- me dical practitioner, and surgeons in particular, with ability 

phy (CT) x-ray systems or by magnetic resonance (MR) l0 plan detailed and complicated surgical procedures using 

imaging systems, either of which is capable of generating totally noninvasive diagnostic methods. The dramatic 

three dimensional arrays of data representative of one or images generated by the invention show every evidence of 

more physical properties at various locations within a three being ^ great an improvement in the medical imaging arts 

dimensional volume. More particularly, the '876 patent is as compil ted axial tomography and nuclear magnetic reso- 

directed to a system and method for the display of medical nance imaging. 

images so as to obtain images and representations of internal While the sy&tem and method of the , g76 patent find& its 

bodily structures. The images generated m the practice of the greatcst utilization in analysis and display of CT and MR 

'876 patent provide three-dimensional data for examination imaging daU) the presenl invention is cqiially applicable to 

by physicians, radiologists, and other medical practitioners. systeras employing ultrasound, positron emission 

In conventional x-ray systems, a two-dimensional shadow tomography, ECT (emission computed tomography) and 

image is created based upon the different x-ray absorption 35 MMI (multimodality imaging). Moreover, while particularly 

characteristics of bone and soft tissues. A great improvement applicable to the construction of medical images, the inven- 

on the conventional x-ray system as a diagnostic tool has t i on ^ a is 0 applicable to the display of interior three dimen- 

been provided by CT systems which are x-ray based and s i ona ] sur f aC e structures for any system capable of gener- 

initially were used to produce single two dimensional views at i ng tare e dimensional data arrays in which signal patterns 

depicting transverse slices of a body, object, or patient being 4Q representing the value of at least one physical property 

examined. Three dimensional information was thereafter associated with points in a solid body are present, 

gleaned from CT scan data by generating data for a number A part i cu i ar advantage of the invention is its capability to 

of contiguous slices and using the inferential abilities of the prov id e the medical practitioner with ability to interact with 

radiologist to suggest a three dimensional representation for the system appar atus in real time. Systems which do not 

the various internal organs. In one embodiment of the '876 45 permit interactive use suffer a significant disadvantage since 

patent, shaded and contoured three-dimensional images are a real time display methodology is required for optimal 

generated from the three-dimensional array of data gener- human interaction with the system, particularly for a surgeon 

ated by a sequence of such contiguous CT scans. In the same planning a difficult procedure. For example, in transplant 
way, the newer MR imaging technology is also capable of surgery, it is often difficult to ascertain beforehand the 

generating three-dimensional arrays of data representing 50 pre cise shape or size of a body cavity which is to receive the 

physical properties of interior bodily organs. implant. This is true whether or not the implant comprises 

Moreover, MR systems have the capability to better human tissue or a mechanical device. It therefore would be 

discriminate between various tissue types, not just between important for a surgeon to be able to display the cavity in 

bone and soft tissue. MR imaging systems are also capable three dimensional form and to rotate and section the cavity 

of generating physiological data as well as image data. 5S image at will, before undertaking any invasive procedure. It 

However, whether MR or CT systems are employed, the is also important that the images generated be sharp and 

data have been made available only as a sequence of slices, exhibit excellent contrast. The images generated should also 

and systems have not generally been available for providing depict surface texture wherever possible, 

true three-dimensional images. Display of three-dimensional graphic images on a cathode 

The '876 patent describes three-dimensional data gener- 60 ray tube (CRT) screen has principally been driven by the 

ated either by a CT scanning system oj^byan MR imaging goals and directions of computer aided design (CAD) and 

system, which may b e d isplayed and mOy^^ _ computer aided manufacture (CAM). Systems have been 

. of ways so as to produce^m"g~ ms^la'v screen_or otherd evice developed for displaying solid bodies and for manipulating 

jTmultitude of anaToTntCal teatures that are selectai ik^th e images in desirable fashion to create solid models for 

Viewer's choice. I n a preferred embodiment of the invention, 65 manufactured parts and for rotating and viewing these parts 

the data used to produce the three dimensional images are from a multiplicity of directions. In particular, CAD/CAM 

typically acquired once and then used and re-used to gen- systems have been developed which accept so called wire- 
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frame data. In a wire-frame display format, the display polyline of shorter length. This process of analysis, trian- 

processor is provided with a sequence or list of three- gulation and modification is iterative and is continued until 

dimensional points representative of the end points of line a shortest possible polyline is found between the first and 

segments. These line segments are joined to represent vari- second points. This shortest polyline may be used, for 

ous surfaces. An advantage of these wire frame images is the 5 example, to define a planar cut of the surface going through 

ability to rapidly rotate the image about various axes to the two points between which the shortest path has been 

obtain different views. This is a computationally burden- determined. 

some task which has recently been alleviated through utili- Finding the shortest path on a non-flat surface is a 

zation of high speed, large-scale integrated circuit devices. requirement encountered in numerous practical applications, 

A wire-frame image, even one which has been processed io ei g t) building the shortest road connecting two locations in 

to eliminate hidden lines, may typically comprise a list of a mountainous region, or flying the shortest route between 

50,000 vectors which is displayed on a screen, each "vector*' two locations in a mountainous region while maintaining 

being a (directed) line segment drawn between two points on constant elevation above the ground (typical for helicopter 

a CRT form of display device. More sophisticated graphics military missions). In the first example there may be other 

processing systems not only accept wire-frame data, but also 15 criteria apart from the road length, such as maximum incline 

perform functions such as hidden line removal and contimi- of the road and water crossings that may require some 

ous shading in various colors and/or shades of gray. In such modification to the presented method while, in the second 

systems, the viewing angle is selected by the operator. In example the method is applicable without any modification, 

systems displaying shaded images, the normal vector to the In both examples, solution is possible due to the availability 

surface being displayed varies from point to point on the 20 of digital elevation maps that provide polygonal represen- 

surface and the shading is determined from this normal tation of the earth surface. 

vector and the viewing angle. Thus the information provided Another application is to find the minimal distance 

by the normal vector is critically important in applying the between two points on the surface of the brain (the cortex) 

shading to what is in actuality a two dimensional CRT screen as extracted from an MR data set. This distance is of use to 

image. It is the shading of this image (based originally on 25 reS earchers studying physiological information obtained, for 

wire-frame data) that creates the illusion of three dimen- instance, from functional MR data. The pure straight-line 

sions. Providing these three dimensional clues (shading) to distance is not as useful as the shortest path along the surface 

the human visual system is achieved particularly well by the s j nce tne distance along the surface can be related to 

present invention. The shading data are generally produced physiological activity in the brain, 

from wire-frame data in a graphics processing system which 3° 

converts the vector format of wireframe images to the BRIEF DESCRIPTION OF THE DRAWINGS 

so-called raster scan format that is commonly used in nG x fc a d[ illustrating the three adjacent points 

displaying television images. of a polyline; 

Accordingly, there is a need for improving the computa- nG % ^ & m iUustrati a [n[ intcrsection 

turn and display of three-dimensional informaUon. One of & ^ ^ a common and 

the computational goals is to use a method for finding the _ „ . .„ . , , 

shortest path between two points on a polygonal surface. * 3 15 a dia S rara ^traUng the * ho ^ P 0 ^ 1 ™- 

Such method would be useful in conjunction with CT DETAILED DESCRIPTION OF INVENTION 
scanners, ultrasound devices, MR imaging systems, and 

other systems capable of generating three dimensional data 40 Given two P oints on a surface, it is necessary to define a 

representative of one or more physical properties within a polyline P lying on the surface and passing through these 

body to be studied. These techniques also have application two P oims and > for example, defining a planar cut of the 

in graphics systems for medical images, which are capable surface passing through these two points, 

of interactive use and, at the same time, produce high quality FIG. 1 illustrates a polyline lying on a polygonal surface 

images that provide textural, shading, and other visual clues 45 10 and passing through the two points. Thus, (p a , p 2 . . . ,p nl ) 

to the user. Such graphics systems provide medical practi- is a polyline defined as P^ with vertices (p 1( p 2 . . . ,p Nl ) all 

tioners with the ability to emulate surgical procedures lying on the surface and, with the possible exception of the 

graphically prior to undertaking invasive measures. first and last points, also lying on the edges of a polygonal 

cTTxynvjAnv^cTiTc fvn/EwrmM mesn defining the surface. It is assumed, without loss of 

SUMMARY OF THE INVENTION 50 gcncralityf th * al i polygonal facesof the me sh are triangles. 

In a preferred embodiment of the invention, the shortest if the first and last points of the polyline are not on the edges 

path is found between first and second points through which 0 f the polygonal mesh, the faces on which they lie can be 

is passed a polyline lying on a polygonal surface. The triangulated so that the start and the end points of the 

polyline is a one-dimensional cell consisting of one or more polyline become vertices of the polygonal mesh. L(P) is 

connected lines, and is defined by an ordered list of n+1 55 denoted the length of polyline P. 

points, where n is the number of lines in the polyline and Designating the first point of the polyline as A and the last 

each pair of points (i, n-1) defines a line. The polyline is poim of the polyline ^ B) the polyline is modified such that 

analyzed on a polygonal mesh defining the surface to it win pass through points A and B and have the locally 

determine points lying on the polyline and on edges of the shortest length (become a geodesic on the surface). This 

mesh. The polygonal faces of the mesh are assumed, without 60 mo difi ca tion is an iterative process that starts with the first 

any loss of generality, to be triangles. If the start and end three points of the polyline and creates a new polyline with 

points of the polyline are not on the edges of the mesh, the a shoftcr lcngth that has point pj replaccd by a 0 f points 

faces of the polygonal surface on which the start and end ^ ^ ^ ^ process ^ described in detail below. 

points of the polyline lie are triangulated so that the start and ™ , r . n , 2 2 2\ t u 

... ^ r.i 1 1 1 r™ The new polyline is P-e(p/ p, 1 ' p/) meeting the 

end points become vertices of the polygonal mesh. The 65 con£ jj^ ons . 

polyline is then modified such that it will pass through the 

first and second points of the polyline, creating a new A«p x 2 , B~p n2 2 > L{Pi)<L{P^ 
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The process is repeated for points p*, p 2 2 , p 3 2 , resulting support library or in any computational tool evaluating 

in a new polyline P 3 -(pj 3 , p 2 3 , ■ • ■ »P„3 3 ) meeting the surfaces within a three-dimensional display system. Such 

condition: computational tools may be included in medical modalities 

such as magnetic resonance, CT or X-ray equipment. In 

Atp^, B»p n ^ t L(Pj)<L(P^). s addition, the method may be used in situations where a 

polygonal surface is used as a polygonal representation of 

Eventually, the situation is reached in which the second industrial CAD systems, 

point in the polyline will not move/split. The process pro- While the invention can be embodied in the form of 

ceeds to the second, third and fourth points of the extant computer-implemented processes and apparatuses for prac- 

polyline and repeats the optimization until the second point ^ ticing those processes, it can also be embodied in the form 

(the third if counting from the beginning) stops moving/ of computer program code containing instructions recorded 

splitting. These iterations are continued to the last triple of in tangible media, such as floppy diskettes, CD-ROMs, hard 

points. After that the process starts again with the first triple drives, or any other computer-readable storage medium, 

of points, etc. The iteration process generates a sequence of wherein, when the computer program code is loaded into 

polylines Pj,P 2 ,... that start and end at A and Band satisfy and executed by a computer, the computer becomes an 

the conditions: apparatus for practicing the invention. When implemented 

on a general-purpose microprocessor, the computer program 

i(^i)>i(^2)>^(^3)> code segments configure the microprocessor to create spe- 
cific logic circuits. 

Reprocess stops when it is not possible to move/split the whik Qnly preferrcd features of the invention 

second point after finishing the last triple. At this time the 20 haye been illustrated and describ ed, many modifications and 

process has found the polyline with the locally shortest changes win Qccur tQ those [n mc art }{ ^ thcrcfore> 

length that connect points A and B and lies on the surface. tQ be understood that the appende d claims are intended to 

The move/split procedure can be understood by consid- covef ^ such modifications and changes ^ fall within the 

enng three adjacent points of the polyline Pl , p 2 , p 3 shown tnje - it of the invention> 

in FIG. 1. Point pj=A is a vertex of triangle t^pjNM and 25 W hat is claimed is- 

point p 2 is on the edge of triangle t^-NRM While points A i. a method for finding the shortest path between two 

and B are at vertices of the triangles, all other points are on ims Qn a polygonal surface comprising: 

the edges/vertices of the triangles. Triangle t, is rotated . , , r , . . , - , 

°. ' 4 . . . f, , ■ i * -m. establishing a polyline lying on said surface and passing 

around its edge to align it with the plane of triangle t,. The ,, u c. * - * a a • ♦ a ^ 

, t , i i , b . 4 *\ . . ■ i . i ■ ,„ through a first point and second point on said surface; 

shortest path between points p, and p 3 is a straight line 30 , i LJ ^. 

connecting these two points in the plane! and q is a point of sai ^ P ol y hne on ? ?° meS * de ?™S Sai * 

intersection of this line and the common edge. Then Uv» q. l ° P° ints ^ °° said and 

P,)<L(P., P 2 , P 3 ). If q ^ inside the edge, the process stops. ? n ^ges of said mesh, said mesh having polygonal 

Otherwise, it is necessary to consider what happens when fac f 115 the form of tnan § les; 

point q moves to the end vertex of the edge and beyond, as 35 modifying said polyline such to pass through said first and 

shown in FIG 2 second points of said polyline and thereby create a new 

' For q at the end point of the edge the inequality still holds, P ol y line of shorter len 6 th than said first polyline; and 

but the process determines if the path can be shortened even repeating the analysis, tnangulation and modification 

more by moving the point beyond the vertex. The planar iteratively until a shortest possible polyline is found 

angle Zp^ is less than Jt, otherwise point q would not be *o between said first and second points, 

the end point of the edge common to triangles t^qM and 2 - The method of claim 1 wherein start and end points of 

t 2 =qRM. If all triangles of the mesh sharing the vertex q and said polyline lie on other than an edge of said polygonal 

connecting triangles l x and t^ are flattened on the plane mesh and including, prior to modifying said polyline: 

together with t n and a new angle Z.pjqp 3 is obtained. The triangulating faces of said polygonal surface upon which 

angle value may be different though the same lines p 3 q and 45 said start and end points of said polyline lie to become 

qp 3 are used for angle definition. The lines may move after vertices of said polygonal mesh, 

flattening all triangles around vertex q. If this angle is greater 3. The method of claim 1 wherein one of said start and end 

then 7i the process stops here. Otherwise the process can split points of said polyline lies on other than an edge of said 

point q into several points, one for each of the edges of the polygonal mesh and including, prior to modifying said 

triangles that share q. The process results in a polyline p 1? 50 polyline: 

<h» • ■ • Mk> p3 while, if there is just one triangle between tj triangulating a face of said polygonal surface upon which 

and t2, and p 3 is a vertex, then k may be 0. The shortest said one of the start and end points of said polyline lies 

polyline that connects points p a and p 3 and lies inside all to become a vertex of said polygonal mesh, 

these triangles is part of a convex hull of points p v Q u . . . , 4. The method of claim 1 wherein said polygonal surface 

Q^ps, where points Q a , . . . y Q k are ends of edges emanating 55 is a polygonal representation of one of the group consisting 

from vertex q as shown in FIG. 3. of earth's surface, a computer graphics display, and indus- 

The polyline obtained through this process satisfies the trial CAD parts, 

condition L(p x , q lt . . . ,q*, p 3 )<L(pi, p 2 > P3) and the process 5. The method of claim 1 wherein said polygonal surface 

creates a polyline P^pj 2 , . . . ,p M2 2 )=(Pi> Qi* • ■ • is a polygonal representation of a human body part extracted 

P?> • ■ • ypni) f° r which the condition L(P^}<L (Pj) holds. 60 from a medical modality data set from one of the group 

Because this process makes use of polylines, the memory consisting of magnetic resonance, computerized 

capacity required for storing the points through which the tomography, and X-ray. 

polyline passes is greatly reduced from that which would be 6. A storage medium encoded with a machine-readable 

needed if polylines were not employed. As a result, the computer program for finding the shortest path between two 

process runs faster and the costs associated with implement- 65 points on a polygonal surface, the storage medium including 

ing the process are reduced. The process may be imple- instructions for causing a computer to implement a method 

mented in a computer graphics system including a graphics comprising: 
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establishing a polyline lying on said surface and passing 
through a first point and a second point on said surface; 

analyzing said polyline on a polygonal mesh defining said 
surface to determine points lying on said polyline and 
on edges of said mesh, said mesh having polygonal s 
faces in the form of triangles; 

modifying said polyline to pass through said first and 
second points of said polyline and thereby create a new 
polyline of shorter length than said first polyline; and 

repeating the analysis, triangulation and modification 
iteratively until a shortest possible polyline is found 
between said first and second points. 

7. The storage medium of claim 6 wherein start and end 
points of said polyline lie on other than an edge of said 
polygonal mesh and including, prior to modifying said 
polyline: 

triangulating faces of said polygonal surface upon which 
said start and end points of said polyline lie to become 
vertices of said polygonal mesh. 20 

8. The storage medium of claim 6 wherein one of said 
start and end points of said polyline lies on other than an 
edge of said polygonal mesh and including, prior to modi- 
fying said polyline: 

triangulating a face of said polygonal surface upon which 25 
said one of the start and end points of said polyline lies 
to become a vertex of said polygonal mesh. 

9. The storage medium of claim 4 wherein the polygonal 
surface is a polygonal representation of one of the group 
consisting of earth's surface, a computer graphics display, 30 
and industrial CAD parts. 

10. The storage medium of claim 6 wherein said polygo- 
nal surface is a polygonal representation of a human body 
part extracted from a medical modality data set from one of 
the group consisting of magnetic resonance, computerized 35 
tomography, and X-ray. 

11 . A computer program for use with a graphics display 
device, said computer program comprising: 

a computer usable medium capable of having computer 
readable program code means embodied in said 40 
medium for finding the shortest path between two 
points on a polygonal surface; and 

computer readable program code embedded in said 
medium and adapted to operate a computer to perform 
the steps of: 

establishing a polyline lying on said surface and pass- 
ing through a first point and a second point on said 
surface; 

establishing a polyline lying on said surface and pass- 
ing through a first point and a second point on said 
surface; 

analyzing said polyline on a polygonal mesh defining 
said surface to determine points lying on said 
polyline and on edges of said mesh, said mesh 
having polygonal faces in the form of triangles; 55 

modifying said polyline to pass through said first and 
second points of said polyline and thereby create a 
new polyline of shorter length than said first 
polyline; and 

repeating the analysis, triangulation and modification 60 
iteratively until a shortest possible polyline is found 
between said first and second points. 

12. The computer program of claim 11 wherein start and 
end points of said polyline lie on other than an edge of said 
polygonal mesh and including, prior to the step of modifying 
said polyline: 
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triangulating faces of said polygonal surface upon which 
said start and end points of said polyline lie to become 
vertices of said polygonal mesh. 

13. The computer program of claim 11 wherein one of 
said start and end points of said polyline lies on other than 
on edge of said polygonal mesh and including, prior to the 
step of modifying said polyline: 

triangulating a face of said polygonal surface upon which 
said one of the start and end points of said polyline lies 
to become a vertex of said polygonal mesh. 

14. The computer program of claim 11 wherein the 
polygonal surface is a polygonal representation of one of the 
group consisting of earth's surface, an image upon said 
computer graphics display, and industrial CAD parts. 

15. The computer program of claim 11 wherein said 
polygonal surface is a polygonal representation of a human 
body part extracted from a medical modality data set from 
one of the group consisting of magnetic resonance, comput- 
erized tomography, and X-ray. 

16. An article of manufacture for use in a computer having 
an operating system and a graphics support library for 
finding, or a graphics display, the shortest path between two 
points on a polygonal surface, 

said article of manufacture comprising a storage medium 
having computer readable program code embedded 
therein, said program code being adapted to operate 
said computer to perform the steps of: 
establishing a polyline lying on said surface and pass- 
ing through a first point and a second point on said 
surface; 

analyzing said polyline on a polygonal mesh defining 
said surface to determine points lying on said 
polyline and edges of said mesh, said mesh having 
polygonal faces in the form of triangles; 

modifying said polyline to pass through said first and 
second points of said polyline and thereby create a 
new polyline of shorter length than said first 
polyline; and 

repeating the analysis, triangulation and modification 
iteratively until a shortest possible polyline is found 
between said first and second points. 

17. The article of manufacture of claim 16 wherein start 
and end points of said polyline lie on other than an edge of 
said polygonal mesh and including, prior to the step of 
modifying said polyline: 

triangulating faces of said polygonal surface upon which 
start and end points of said polyline lie to become 
vertices of said polygonal mesh. 

18. The article of manufacture of claim 16 wherein said 
one of said start and end points of said polyline lies on other 
than an edge of said polygonal mesh and including, prior to 
the step of modifying said polyline: 

triangulating a face of said polygonal surface upon which 
said one of the start and end points of said polyline lies 
to become a vertex of said polygonal mesh. 

19. The article of manufacture of claim 16 wherein the 
polygonal surface is a polygonal representation of one of the 
group consisting of earth's surface, said computer display, 
and industrial CAD parts. 

20. The article of manufacture of claim 16 wherein said 
polygonal surface is a polygonal representation of a human 
body part extracted from a medical modality data set from 
one of the group consisting of magnetic resonance, comput- 
erized tomography, and X-ray. 
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